Preprint typeset in JHEP style - HYPER VERSION 



The conformal anomaly of k-strings 



Pietro Giudice, Ferdinando Gliozzi, Stefano Lottini 

Dipartimento di Fisica Teorica, Universita di Torino and 

INFN, Sezione di Torino 

via P.Giuria 1, 1-10125 Torino, Italy 

giudice , gliozzi , lottini@to . inf n . it 



Abstract: Simple scaling properties of correlation functions of a confining gauge theory 
in d-dimensions lead to the conclusion that k-string dynamics is described, in the infrared 
limit, by a two-dimensional conformal field theory with conformal anomaly c = (d — 2)ak/cr, 
where is the k-string tension and a that of the fundamental representation. This result 
applies to any gauge theory with stable k-strings. We check it in a 3D Z4 gauge model at 
finite temperature, where a string effect directly related to c can be clearly identified. 
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1. Introduction 

One of the simplest and most general consequences of the effective string description of the 
quark confinement [|l|, ^, [| is the presence of measurable effects on physical observables of 
the gauge theory, produced by the quantum fluctuations of that string ||, ||. The most 
widely known is the Liischer correction to the interquark confining potential V at large 
distance r 

V(r) = ar + 2fi-(d- 2)^- + 0(l/r 2 ) , (1.1) 

where a is the string tension and fi a self-energy term. The attractive, Coulomb-like 
correction is universal in the sense that it is the same for whatever confining gauge theory 
and depends only on the d — 2 transverse oscillation modes of the string. 

A similar universal effect has been found in the low temperature behaviour of the string 
tension || 

a(T) = 0- (d-2)-T 2 + 0(T 4 ). (1.2) 
6 

Both effects may be rephrased by saying that the infrared limit of the effective string 
is described by a two-dimensional conformal field theory (CFT) with conformal anomaly 
c = d — 2. In this language, the (generalised) Liischer term — is the zero-point energy 
of a 2D system of size r with Dirichlet boundary conditions, while the — ^-T 2 term is 
the zero-point energy density in a cylinder (i.e. the string world-sheet of the Polyakov 
correlator) of period L = l/T |7|. 

In this paper we propose a method to extend these results to a more general class of 
confining objects of SU(iV) gauge theories, the k-strings, describing the infrared behaviour 
of the flux tube joining sources in representations with N— ality k. There are of course 
infinitely many irreducible representations corresponding to the same value of k. No matter 
what representation St is chosen, the stable string tension cr& depends only on the iV— ality k 
of i.e. on the number (modulo N) of copies of the fundamental representation needed to 
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build St by tensor product, since the sources may be screened by gluons. As a consequence, 
at sufficiently large distances, the heavier strings find it energetically favourable to decay 
into the string of smallest string tension, called k-string. 

The spectrum of k-string tensions has been extensively studied in recent years, in 
the continuum ]8|-[|l6[] as well as on the lattice [17|-|24j]. So far, in numerical analyses 



one typically measured the temperature-dependent k-string tensions (Jk(T) through the 
Polyakov correlators and then extrapolated to T = using ( |1.2[ ), hence assuming a free 
bosonic string behaviour. 

Recently this assumption has been questioned by a numerical experiment. It showed 
that in a 3D Z4 gauge theory, though the 1-string fitted perfectly the free string formulae 
with a much higher precision than in the SU(A r ) case, the 2-string failed to meet free string 



expectations [25|. One could object that there is no compelling reason for a 2-string of a 
Z4 gauge system to behave like a 2-string of SU(iV) gauge system; a k-string can be seen 
as a bound state of k 1-strings and the binding force would presumably depend on the 
specific properties of the gauge system. 

On the other hand, from a theoretical point of view there are good reasons to expect 
values of c larger than d— 2. In fact the conformal anomaly can be thought of as counting the 
number of degrees of freedom of the k-string. Therefore the relevant degrees of freedom 
are not only the transverse displacements but also the splitting of the k-string into its 
constituent strings. If the mutual interactions where negligible, each constituent string 
could vibrate independently so we had c = k (d — 2). Thus we expect that c can vary in 
the range 

d-2< c <k(d-2). (1.3) 

An unexpected insight into the actual value of c comes when considering the infrared 
properties of the N-point Polyakov correlators related to the baryon vertex of SU(A^) 

(P f (ri)P f (r 2 )...P f (f N )) T , (1.4) 

where Pf{fi) is the Polyakov line in the fundamental representation identified by d — 
1 spatial coordinates r*j and directed along the periodic temporal direction of size L = 
1/T. If all the distances |r*j — fj\ are much larger than any other relevant scale, this 
correlator is expected to obey a simple scaling property. When combining this fact with 
the circumstance that, depending on the positions fi of the sources, some strings of the 
baryon vertex may coalesce into k-strings [^, [l3| , one obtains the geometric constraint 

a k (T) = a k -(d-2)^T 2 + 0(T 3 ) , (1.5) 

which is the main result of this paper. We check this formula in a 3D Z4 gauge model 
where, thanks to duality, very efficient simulation techniques are available, yielding high 
precision results which give fairly convincing evidence for the scaling law ( |1.5| ). 

The contents of this paper are as follows. In the next Section we expose in detail 
the above-mentioned scaling argument, while in the following Section we describe a lattice 
calculation on a three-dimensional Z4 gauge theory where combining duality transformation 
with highly efficient simulation techniques it is possible to confirm that the stable 2-string 
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matches nicely Eq. (|1.5| ), We finish with a discussion of our results and some of their 
implications. 



2. Scaling form of the Polyakov correlators 



The main role of the string picture of confinement is to fix the functional form of the 
vacuum expectation value of gauge invariant operators in the infrared limit. It predicts 
two different asymptotic behaviours of the correlation function of a pair of Polyakov loops 

ri| and L = 1/T are large, depending on the value 



(Pf(ri)Pj(r2)T> when both r = \ f\ 



of the ratio 2r/L. Using (1.1) and (1.2) we can write 



- z log(P / (r 1 )P / T (r - 2 ) )T 



V(r) + 0{e- 7vL/r ) , 2r < L , 



L 



log(P f (n) P}(r 2 )) T 



1 9r 

a(T)r + 2^+ — log j- + 0(e^ r / L ) , 



2r > L . 



(2.1) 



(2.2) 



where V(r) and cr(T) are given by (U) and (|T|). There are strong indications that the 
r~ 3 and T 3 terms are zero and the r 4 and T are universal (see for instance the discussion 
in and references quoted therein). For our purposes we need only the first universal 
term, which is directly related to the central charge of the CFT describing the IR limit 
of the underlying confining string. In this approximation the Polyakov loop correlation 
functions should decay at large r while keeping constant T as 



(P f (n) Pj(f 2 ))T « exp [— <r(T) r L — 2fiL] . 



(2.3) 



Similar expansions are expected to be valid also for Polyakov correlators describing more 
specific features of SU(iV) gauge theory, like those involving baryonic vertices. 

A baryon vertex is a gauge-invariant 
coupling of iV multiplets in the fundamen- 
tal representation / which gives rise to 
configurations of finite energy, or bary- 
onic potential, with ./V external sources. 
At finite temperature T these sources are 
the Polyakov lines Pt(fi). Assume for a 
moment N = 3. When the mutual dis- 
tances \fi — fj\ = 1,2,3) are all large, 
three confining strings of chromo-electric 
flux form, which, starting from the three 
sources, meet a common junction; their 
world-sheet forms a three-bladed surface 
with a common intersection [28, ^] (see 

Fig-D- 




Figure 1: The three-bladed world-sheet of a static 
SU(3) baryon. 



The total area A of this world-sheet is L£(f\, f^, 7*3) where £ is the total length of the 
string. The stable configuration (hence the position of the common junction) is the one 
minimising A. The balance of tensions implies angles of 4^ between the blades. 
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The complete functional form of the 3-point correlator is substantially unknown. 
Nonetheless, the static baryon potential, defined as 



V = — lim T log( P f (n) P f (f 2 ) P f (f 3 ) ) T , 

X — >oo 



(2.4) 



has a simple form in the IR limit: 



V = cr£(ri,f 2 ,r 3 ) + 3/i+ ... 



(2.5) 



The universal 1/r corrections have been calculated in 

In the IR limit at finite temperature, i.e. — rA ^> L\/i ^ j, we assume, in analogy 
with dpi), 



(P f (n) P f {r 2 ) . ..P f (f N ) ) T = e~ FN ~ exp [-cr(T) L£(fi, r 2 , . . . , r N ) - Nfi L] , (2.6) 



or, more explicitly, 

F N (i,L) 



oLi-{d- 2) J- + iV/iL , (\fi - fj\ »LV^ j) 
b L 



(2.7) 



where the coefficient of the l/L term specifies that in this IR limit the behaviour of the 
baryon flux distribution is described by a CFT with conformal anomaly d — 2 on the string 
world-sheet singled out by the position of the external sources. 

The above considerations are completely standard and nothing new happens so far. 
The surprise comes when considering the 

latter expression in the case N > 3. No- ' / > 

tice that, depending on the positions r; L of / / 

the sources, some fundamental strings of the ^, „---'" r 

baryon vertex may coalesce into k-strings 
|26| . As a consequence, the shape of the 'v- 
world-sheet changes in order to balance the 
string tensions and £ becomes a weighted 
sum, where a k-string of length A contributes 
with a term A , where cr(T) is given by 



1.2(), while 

<7fc(T) = 



Ok ~ c k \T 2 + 0(T 3 ) , 
b 



where is the conformal anomaly of the 
k-string. When the string is in the funda- 
mental representation the T 3 term (fT^) is 
missing [pi], 



Figure 2: Confining string configuration of a 
SU(4) static baryon. The thick line is a 2-string. 



tension ratio is 



33], however we cannot exclude it in the k-string with k > 1. The string 



a a 



d-2 



2L)1 T 2 + 0(T 3 
J 6 



a(T) 

if the coefficient of the T 2 term does not vanish the asymptotic functional form (^] 
the free energy gets modified. 



(2- 



of 
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As a simple, illustrative example, let us consider the chromo-electric flux distribution 
of a 4D SU(4) gauge system generated by four external quarks placed at the vertices of 
a regular tetrahedron (see Fig. ^). Preparing the external sources in a symmetric con- 
figuration does not necessarily imply that the distribution of the gauge flux preserves the 
tetrahedral symmetry. In fact, the formation of a 2-string breaks this symmetry (see Fig.Q). 
The actual symmetry breaking or restoration depends on the cost in free energy of the con- 
figuration, of course. It is easy to show that, when < -^=, the tetrahedral symmetry is 
spontaneously broken [13] and the length A of the 2-string, which is a function of the ratio 
(72 /cr, is given by 



A 



r (72 



y/2 -a- 2 

where r is the edge length. The free energy F4 of the four-quark system is 



(2.10) 



F 4 = cr(T) L £{T) + 4/xL = a(T) L I 2r\ 



\2a(T)J v/2 a(T) ' 



Now the total length £ of the string may depend on T through the ratio a 2 {T) / a{T). 
Expanding in T = 1/L as in ( |2.9[ ) we get 

,*(0) _ a (d-2,_C2_\ A" 
a 2 J L 



F 4 = aL£(0)-^ 
b 



(d-2) 



L 



a 



+ 4fiL + 0(l/L 2 ) 



(2.12) 



Clearly the term proportional to A violates the expected asymptotic form of the free energy 
(pl|), unless c 2 = {d- 2) <7 2 /cj . 

More generally, the baryonic free energy 
keeps the expected asymptotic form ( |2.7| ) only 
if the world-sheet shape does not change while 
varying T. Now in a generic string configura- 
tion contributing to SU(iV) baryon potential, 
the angles at a junction of three arbitrary k- 
strings are given by (see Fig.||) 

a^Tf + auiTf-a^Tf 




cos 8j 



2a j (T)a k (T) 



(2.13) 



Figure 3: The balance of the string tensions. 



and others obtained by cyclic permutations of 

the indices k. As a consequence, these angles are kept fixed only if all the string tension 
ratios are constant up to T 3 terms, i.e. only if 

Ci_ _ Cj_ _ Cfc _ (d - 2) 

(7,; (7i <7fc 



(7; 



(7 



(2.14) 



which leads directly to 



<7k(T) = Ok - (d- 2) 



^T 2 + 0(T 3 ) 

6 (7 



(2.15) 



as anticipated in the Introduction. 



- 5 - 



3. The 3D Z 4 gauge model and its dual 



The above general argument on the finite temperature corrections of the k-string tensions 
suggests a different behaviour with respect to the usual assumption that these corrections 
are those produced by a free bosonic string. Since the comparison with theoretical predic- 
tions of k-string tensions is sensitive to this behaviour, it is important to check its validity. 

In this Section we address such a question with a lattice calculation in a 3D Z4 gauge 
theory which is perhaps the simplest gauge system where there is a stable 2-string. 

We work on a periodic cubic lattice L x L x L T . The degrees of freedom are the fourth 
roots of the identity Q (£ 4 = 1), defined on the links I of the lattice. The partition function 
is 

Z(Pf,Pff) = H E (3.1) 
I Q=±i,±i 

where the sum is extended to all plaquettes p of the lattice and U p = ILepCz; Pf anci Pff 
are two coupling constants. When they vary in a suitable range the system belongs to a 
confining phase. In analogy with the SU(iV) case we say that U p is in the fundamental 
(/) representation while U p lies in the double-fundamental (//) representation. From a 
computational point of view it is convenient to recast Z as the partition function of two 
coupled Z2 gauge systems 

m>Pff) = U E e £^ + ^ + /W P , ( U p = 11^, V p = Hv] . (3.2) 

I {Ui=±l,V t =±l} \ lap iGp J 

The external sources generating the 1-string and the 2-string are given respectively by the 
two products 

P f (r) = J] Ui or H V t ; P ff (r) = J] , (3.3) 

where 7^ is a closed path in the lattice which winds once around the temporal direction L T 
and passes through r. 

This model, as any three-dimensional abelian gauge system, admits a spin model as 
its dual. We have recently shown that this Z4 gauge model is dual to a spin model 
with global Z4 symmetry which can be written as a symmetric Ashkin- Teller (AT) model, 
i.e. two coupled Ising models defined by the two-parameter action 

Sat = - E [Pi a ^ a y + T ^ T y) + ot(<x x a y r x T y )] , (3.4) 

(xy) 

where a x and t x are Ising variables (a, r = ±1) associated to the site x and the sum is over 
all the links (xy) of the dual cubic lattice. The two couplings a and /? are related to the 



gauge couplings j3f and (3ff by |34] 



1 / (coth/?j + tanh/3j tanh/3jj)(coth /3j + tanh/3j coth/3jj) \ 
a= 4 l0g V 2 + tanh/3 // + coth/3 // ~)> (3 ' 5) 



1 fl + tanh 2 8 f tanh (3 f f 

— lOg X " 

4 V tanh 2 /?/ + tanh% 
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The duality transformation maps any physical observable of the gauge theory into a cor- 
responding observable of the spin model. In particular it is well known that the Wilson 
loops are related to suitable flips of the couplings of the spin model. We have found the 
identity 

( V P ) gauge = ( e-W+^ry)*^ (3 7) 

generalising the known dual identity of the Ising model. Similarly, flipping the signs of 
both spins a x and r x we get the plaquette variable in the k = 2 representation as ( UpVp ). 
Combining together a suitable set of plaquettes we may build up any Wilson loop or 
Polyakov-Polyakov correlator with k = 1 or k = 2. 



3.1 Monte Carlo simulations 



Our interest in writing this model in terms of Ising variables is twofold. 

First we can easily apply for the simulation a very efficient non-local algorithm p5| . 
basically similar to the standard Fortuin-Kasteleyn cluster method of Ising systems: each 
update step is composed by an update of the a variables using the current values of the 
r's as a background (thus locally changing the coupling from (3 to (3 ± a according to the 
value of T x T y on the link (xy)), followed by an update of the r's using the a values as 
background. 

Secondly, by nipping a suitable 
set of couplings, we can insert any 
Wilson loop or Polyakov correlator di- 
rectly in the Boltzmann factor, pro- 
ducing results with very high preci- 
sion. If, for instance, we generate a 
sequence of Monte Carlo configura- 
tions where the a x couplings of all 
the links crossing the cylindric sur- 
face bounded by the Polyakov lines 
7^ and 7^ 2 are flipped, then the av- 
erage of whatever observable Q is ac- 
tually the quantity 



{QP f {n)P){r 2 )) T 



(3.8) 
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0.01271(2) 
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Table 1: String tension differences Aa(T) as resulting 
from fits to Eq.( 3.10 ). Aer(0) is evaluated from a fit 



In our numerical experiment we choose to E q-(tLHD usin S AfT ( T ) data - AU the quantities are 

Q = P f {f x ) Pj(f 2 ), therefore the cor- expressed in lattice spacing units. 

responding averages yield directly, ac- 
cording to (|Q|), the ratio 



R(\n-f 2 \,T)=iMmium 

(P f (n)P}m T 



(3.9) 
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We estimated this vacuum expectation value with a very powerful method, based on the 
linking properties of the FK clusters |36|]: for each FK configuration generated by the 
above-mentioned algorithm one looks for paths in the clusters linked with the loops 7^ 
and 7^ 2 . If there is no path of this kind we put Q = 1, otherwise we set Q = 0. The 
algorithm we used to determine the linking properties is described in p7[] . This method 
leads to an estimate of R{\ri — r^l 1 T) with reduced variance with respect to the conventional 
numerical estimates. 



3.2 Results 

We performed our Monte Carlo simulations on the AT model, at two different points of 
the confining region, for which we measured previously the string tensions [^5| (see Table 
1) . We worked on a cubic periodic lattice of size 128 x 128 x N T with N T chosen in such 
a way that temperature of our simulations ranged from T/T c ~ 0.5 to T/T c ~ 0.8 and we 
took the averages over 10 6 configurations in each point. 




20 30 r 40 50 60 

Figure 4: Plot of G(r)ff/G(r)f on a log scale. 



The large distance behaviour of the data is well described by a purely exponential 
behaviour (see FigM) 



(Pff(n)Pff(r 2 )) 



G(r) ff /G(r)f = \UIX^U^L±ILL a e -A«rN T ^ (310) 



T 



with Aff = aff — Of. Comparison with ( p.2| ) shows that the logarithmic term, which 
is a potential source of systematic errors when neglected in Polyakov correlators, here is 



cancelled in the ratio. Since ( 3.10 ) is an asymptotic expression, valid in the IR limit, we 
fitted the data to the exponential by progressively discarding the short distance points and 
taking all the values in the range r m j n < r < r max = 60 a with r m j n varying from 15 to 40 
lattice spacings a. The resulting value of Act turns out to be very stable, as Fig. || shows. 
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Figure 5: The fitted value of Act to Eq.(3.10) as a function of the minimal distance r m i n considered 



for different values of N T . The large plateaux show the stability of the fit. 



The whole set of values of Aa(T) as functions of the inverse temperature N T = 1/T are 
reported in Table 1. 

According to Eq. ( |1.5| ) , in the low temperature limit we expect the asymptotic behaviour 



Act(T) = Act(0) 1 



7T 

6a~ 



T ) + 0(T d 



(3.11) 



Assuming for a the values estimated in [25], we used Acr(0) as the only fitting parameter. 
Neglecting one or two points too close to T c we got very good fits to ( 3,11| ) as shown in 
Figs. ^ and |?]. The fitting parameter A <r(0) as well as the estimates of A a(T) are reported 
in Table 1. Note that Acr(0) agrees with the difference of the previous estimates [25], also 
reported in Table 1, however the error is reduced by a factor of 25 (first set of data) and 
even of 50 (second set of data). The reason of this gain in precision is due to the fact that 
g 2 was evaluated from a fit to ( |2.3|) , even taking in account the Next-to-Leading-Order 
terms, which was rather poor because the 2-string does not behave as a free bosonic string 
[25]. On the contrary our fits to (|3~T0|) and flOl)) are very stable and the corresponding 
reduced x 2 /d-o.f are of the order of 1 or less. 



4. Discussion 



Our calculations in this paper were in two parts. In the first part we investigated the 
string nature of long flux tubes generated in an SU(iV) gauge theory by external sources in 
representations with A— ality k > 1. In a previous work (2f|] we were led to the conclusion 
that these k-strings are not adequately described by the free bosonic string picture. In this 
paper we argued that the central charge of the effective k-string theory is not simply d — 2, 
like in the free bosonic string, but = (d — 2)^-. This simple recipe was a straightforward 
consequence of demanding that the asymptotic functional form of Polyakov correlators 
associated to the baryon vertex should not change when varying the mutual positions of 
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Figure 6: Values of Act(T) from the first set of data of Table f versus T 2 = 1/N?. The solid line 



is a one-parameter fit to Eq.(3.11 
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Figure 7: Same as Fig.0, but for the second set of data of Table f . 



the Polyakov lines, the reason being that certain positions create k-strings inside the baryon 
vertex which modify the functional form of the correlator unless -^jy = constant + 0(T 3 ). 
This in turn fixes unambiguously the value of c\.. Although the geometric derivation is 
quite general and the resulting expression for Ck is appealing for its simplicity, it would be 
very important to find some independent quantum argument in support of it. 

The second part of the paper dealt with 2-strings in a 3D Z4 gauge model and in 
particular with the difference Acx(T) = 02 (T) — <j(T) of the string tensions as a function 
of the temperature. Combining together three essential ingredients, namely the duality 
transformation, an efficient non-local cluster algorithm and finally a choice of flipped links 
which allows to directly measure the ratio of Polyakov correlators belonging to different 
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representations, we obtained values for Ac(T) with unprecedented precision which nicely 



agree with our general formula (1.5). There is however much scope for improving these 
calculations: as Fig. || and [7| show, with little more effort it would be possible to evaluate 
also the corrections of order T 3 and T 4 , that in the case of fundamental string are expected 
to be universal. 
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